Order by disorder and gauge-like degeneracy in quantum pyrochlore antiferromagnet 
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The (three-dimensional) pyrochlore lattice antiferromagnet with Heisenberg spins of large spin length S is 
a highly frustrated model with an macroscopic degeneracy of classical ground states. The zero-point energy 
of (harmonic order) spin wave fluctuations distinguishes a subset of these states. I derive an approximate but 
illuminating effective Hamiltonian, acting within the subspace of Ising spin configurations representing the 
collinear ground states. It consists of products of Ising spins around loops, i.e has the form of a Zi lattice gauge 
theory. The remaining ground state entropy is still infinite but not extensive, being O(L) for system size 0(L 3 ). 
All these ground states have unit cells bigger than those considered previously. 
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What is the nature of the ground state of a highly frustrated 
antiferromagnet with spin length SI For unfrustrated (and 
simply frustrated, e.g. triangular) antiferromagnets, a valid 
recipe is to find the classical ground state (i.e., solution of 
mean-field theory) and dress it with spin-wave fluctuations. 
But many frustrated magnets with vector spins have a large 
classical ground state manifold, having more dimensions than 
the two or three guaranteed by spin-space rotational symme- 
try. So, to address frustrated systems with nontrivial ground 
state degeneracies, one must furthermore identify how quan- 
tum fluctuations elect one of the classical states to be the basis 
of the real ground state (or, alternatively, mix them all, pro- 
ducing a spin-disordered ground state). 

Consider the magnon zero-point energy 1 1 ] at harmonic or- 
der, .F h ({nj}) = J2 m \ hijj m , where oj m runs over all modes 
of spin waves fluctuating around a particular classical ground 
state s,; = S&i with unit vectors {n^}. [Note that !F h is 
defined only on classical ground states: it is just the first 
term, after the classical one, of an expansion in 1/S, and thus 
dominates in the large-5 (semiclassical) limit.] Inequivalent 
ground states do have different spectra {u> m }, so JF h ({tii}) 
does take different values, and singles out a particular ground 
state. In simpler cases with (non massive) degeneracies |2] - 
e.g. face-centered cubic (fee) antiferromagnets - JT h ({ni}) is 
minimized for a unique (modulo rotations) classical configu- 
ration {vii}: this long-range ordered state is our answer (for 
sufficiently large S). Such breaking of a degeneracy by quan- 
tum or thermal fluctuations, or quenched disorder, has been 
called "order by disorder" 0|3|]. 

"Highly frustrated" magnets have similar but larger degen- 
erate manifolds, in that the number of independent angle pa- 
rameters is extensive in the system size. The best-known 
examples are the two-dimensional kagome and the three- 
dimensional pyrochlore lattice, built respectively from corner- 
sharing triangles or tetrahedra. My aim here is to ascertain 
for which classical configuration the quantum fluctuation en- 
ergy T h is smallest; if unique (modulo spin rotation and lat- 
tice symmetries), this configuration would give the long-range 
order of the true ground state. 

This paper will address only the T = ordered state of 
the large-5, pure Heisenberg quantum antiferromagnet with 
nearest-neighbor couplings on the pyrochlore lattice, which 
might describe the low-T state that ZnCr^C^ (S = 3/2) is ap- 
proaching before it undergoes a structural transition |4]. Thus 



our expectation of long-range order does not contradict the ev- 
idence for spin-disordered (spin liquid or valence bond crys- 
tal) states at S = 1/2 |5l|6(|. or the lack of collinear order as 
T -> in the classical model 0H. 

Effective Hamiltonian idea — A major barrier to identify- 
ing the selected state is that there is an infinity of states to 
expand around. This problem is often approached by com- 
puting and comparing the JF h values for special states which 
have exceptional symmetry fioll or a small magnetic cell (IHl . 
For these special ground states, a numerically exact can 
be found by integrating over the magnetic Brillouin zone. Yet 
in principle we need to know JF h for each of the continuously 
infinite ground states, almost all of which are nonperiodic; in- 
deed, the true minimum might not lie in this special subspace 
[In the present case, it has a larger magnetic cell.] 

Instead, my general approach is to derive .F h ({nj}) as 
an effective Hamiltonian | 12], for a generic classical ground 
state, in a crude approximation that has no controlled small 
parameter, but results in an elegant form (products of spins 
around loops). This is faithful in the sense that the approxi- 
mate energies are (mostly) in the same order as the exact ones, 
and it displays clearly which attributes of a configuration af- 
fect its harmonic spin-wave energy. 

The collinear states, in the pyrochlore system, are the spe- 
cial subset of states in which A; = r^n for all % with rji = ±1. 
Apart from a global spin axis, they are simply the ground 
states of the Ising pyrochlore antiferromagnet 11311 : this dis- 
crete subset still has a massive degeneracy, expressed now 
as an extensive ground-state entropy (~ iVln(3/2)). I as- 
sert (based on many checks, but no proof) that every collinear 
state is a local minimum of T h 11411 . in accord with the gen- 
eral behavior in exchange-coupled systems |lL|2|]; I further 
conjecture that the optimal T is attained on a collinear state. 
So we may henceforth limit ourselves to collinear states (an 
enormous reduction of the set) and seek the discrete effective 
Hamiltonian T ({Vi})- 

In the kagome case the discrete coplanar states |9] play 
the same role as collinear states do for the pyrochlore; there, 
it turned out that it every coplanar state has the same spec- 
trum {uj m }, hence J 7 () was independent of the discrete con- 
figuration (among coplanar states). A self-consistent anhar- 
monic calculation was required to reduce this degeneracy to a 
unique state with the final conclusion of long-range order in 
the^Xv^ state (30. Note the independence of the spec- 
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trum depended on the fact that hi ■ hj = —1/2 for neighbor 
spins in every coplanar kagome state. 

On the other hand, in collinear pyrochlore states, hi ■ hj = 
±1, hence the classical degeneracy is split by J rh . The goal 
of this paper is to discover how it is split: what is the (approx- 
imate) analytic form of !F h , what is its energy scale, which 
spin pattern gives the minimum jF h , and how large is the re- 
maining degeneracy. It will be shown that there is an infinite 
(discrete) number of minimum-energy states, but the ground- 
state entropy is not extensive in the system size. 

Hamiltonian and equation of motion — The pyrochlore lat- 
tice is the medial lattice of a diamond lattice, i.e. the py- 
rochlore sites are the midpoints of the diamond bonds. Every 
tetrahedron of sites in the pyrochlore lattice is centered by a 
diamond site. With N unit cells, there are N s — AN spins. 
Our Hamiltonian contains only isotropic, antiferromagnet ex- 
change coupling between nearest neighbor spins s, of length 
S on the pyrochlore lattice: H = J Yltij) s i ' s j> where "{ij)" 
counts each neighbor pair just once. 

Now, define the tetrahedron spin L Q = J^iea s * where a 
(and other Greek indices) runs over the tetrahedron centers, 
and "i £ a" means i is one of the four sites in tetrahedron a. 
Then the Hamiltonian can be rewritten TL = \ Jj2 a an( ^ 
the classical ground states are the (very many) states satisfying 
L a =0. 

The {uj m } which enter !F h are the same as the eigenfre- 
quencies of the (linearized) classical dynamics. The classical 
equations of motion are 

hki = S~ 1 s i x hi = JS~ 1 s i x L Q (1) 

where h, = 5H/Ssi = Jj2(ij) s j = J E a: j £a ( L a ~ s 0- In 
explaining the slow dynamics of the classical system, Ref. [8] 
had noted that Q implies 

ht a = -S- 1 J^2s %{al3) xL f - ! , (2) 



where the sum is over the four tetrahedra (3 that are nearest 
to a, and i(a(3) denotes the (unique) spin shared by a and (3. 
Without loss of generality, take s^ ' = i]iSi in our (assumed) 
collinear state, and parametrize the deviations as <5L Q (t) = 
(SL ax (t),5L ay (t),0): the linearized spin dynamics reduces 
to 

P 

in which a 2N x 2N matrix fi is defined with elements 
fi a/ 3 = r).i( a 0). So, via the trick of using tetrahedron spins, the 
dynamical matrix is the classical spin configuration; Now, if 
we can only massage the formulation so it appears as a pertur- 
bation, the expansion will generate an effective Hamiltonian. 

Trace expansion and loop effective Hamiltonian — Iterat- 
ing l|3} gives the eigenvalue equation 

(^) 2 <5L Q = -h 2 5L a = {SJ) 2 ^(/i 2 ) Q7 <$L 7 . (4) 
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TABLE I: Exact versus approximate zero-point energy, as a multiple 
of JSN S . The magnetic unit cell has -/V mag spins. Last column is 
from Eq. Q with A — 4, from the expansion of {5J to 0(M 4 ). 



State iV ma g (r r ) (r r r s ) 


•F exact A — 4 


a,b 4,8 +1 +1 
c 32 -1 +1 
d 64 


0.4498 0.439 
0.4245 0.414 
0.4460 0.430 



A square root and trace now express the desired effective 
Hamiltonian iffTet ITtIi : 

•^(M) = ij>*» s \ JS Tr (V . (5) 

m 

The matrix has constant diagonal elements, (fi 2 ) aa = 4, so 
let us break it up as fi 2 = 41 + M, where I is the identity ma- 
trix; M has elements of form rjirjj and only connects (nearest) 
even sites of the diamond lattice. Next, insert into (|5jl a for- 
mal Taylor expansion of (41 + M) 1 / 2 Every term Tr(M") is 

a product 11^=1 Vik over a closed walk of 2n steps on the dia- 
mond lattice. A step retraced twice contributes a trivial factor 
r]f = 1. Thus 

r ier 

where T runs over loops (without acute angles) in the py- 
rochlore lattice. Eq. (|5Jl has exactly the form of an Ising lattice 
gauge model fl8ll on a diamond lattice. The shortest loops 
r are 6-step and 8-step loops, which appearing first in the 
Tr(M 3 ) and Tr(M 4 ) terms (respectively) of the expansion. 
Such a hexagon and puckered octagon are outlined in Fig.Q 
call their coefficients and K$ . 

We can rewrite © in a more convenient (for some pur- 
poses) form. Let t t = Jliehexagon r Vi = ±1. *e product 
around one of the hexagons. This is a convenient set of vari- 
ables, since distinct but "gauge" equivalent Ising states map 
to the same configuration of r r 's, and any product around a 
longer loop can be written in terms of r r 's. We place the value 
r r at the center of each its hexagon; these sites form a new py- 
rochlore lattice complementary to the old one (and having the 
same lattice constant). In this notation, becomes fl9ll 

^ h = -^ 6 $> r -iif 8 J> r r s + ... (7) 

r (r,s) 

This has the appearance of an Ising Hamiltonian with an ex- 
ternal field Kq < 0, and a coupling K$ > 0; both terms are 
satisfied by r r = —1, i.e. every hexagon has a negative loop 
product. Configurations of that realize this and satisfy 
the classical ground state condition, J2ie a Vi = - which do 
indeed exist - are the ground states to harmonic order [i.e. to 
O(JS)]. Fig.^c) shows the simplest harmonic ground state 
to construct, with 32 spins in the magnetic unit cell; the small- 
est cell has 16 spins. [All previously studied cases were Q = 
states with a 4-spin magnetic cell, for which the collinear state 
always has Ylier Vi = +1' the worst case.] 
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In Table[I] the numerically exact zero point energy is com- 
pared to taken to order M 4 . The configurations {rji} were 
constructed to have different average values (as listed) of the 
terms in Q. Configurations (a,b,c) are shown in Fig.[T] A fit 
based on just the three exact energies in Table H] would give 
Kf- = -1.3 x 1(T 2 and Kf = +1.8 x 1CT 2 . The 0(M 3 ) 
analytic result is Kf 1 = -3/2 8 = -1.2 x 10~ 2 , and the 
0(M 4 ) term gives Kf = 15/2 11 = +0.7 x 10~ 2 . 

It will be seen that a term in (|6} for a 2n-step loop actually 
gets contributions from infinitely many terms in the expansion 
of (|5}, as terms in M m with m > n can represent the loop T 
plus retraced steps. In a sense, the higher powers of con- 
tain additional contributions of I; one might speculate that a 
different breakup, /x 2 = AAl + M', is more accurate. As 
the spectrum of fi 2 fills the interval [0, 16], one choice is to 
expand around its center, i.e. AA = 8, to optimize the conver- 
gence properties; the original choice AA = 4 is also natural 
since M' has no diagonal terms. As an ad-hoc compromise, 
I suggest the geometric mean A = \f2; exactly that value 
turned out to be optimal in a subsequent, systematic graph- 
ical resummation, using a Bethe lattice approximation Hal . 
[Those numerical energies agreed to with numerically 
exact integration of the zero-point energy for a large sample 
of classical ground states llql .l 

Gauge-like symmetry and ground-state degeneracy — 
Eq. l|5} has an exact "gauge" invariance. Let 77^ ^ = 
UaUfif]i{a^) where u a = ±1 arbitrarily; in matrix notation, 
fx' = U/xU \ where U = diag({u Q }). Thus fi' is similar 
to fi, the two configurations have the same harmonic eigen- 
value spectrum and exactly the same value for T h {{-qi\): in 
other words, this is a generic way to make {77^} degenerate 
with I place the term "gauge" in quotes, since genuinely 
gauge-equivalent configurations are the same state which has 
been labeled redundantly, whereas here they are distinct quan- 
tum states. Note too that we must uphold an independent 
(non-"gauge"-invariant) condition that J^iea Vi = in ev- 
ery tetrahedron, so this spin configuration is still a classical 
ground state. 

We can construct a relatively simple family of ground states 
of JF h . A slice (as marked in Fig.Q has sites with four dif- 
ferent z values, constant along the diagonal lines. Let the rji 
along each diagonal line I be alternating in sign; to fix an 
overall sign for each, choose one (110) and one (110) plane 
normal to the lines, and let vg be the spin in that layer. Then, 
the extreme tip sites (as seen in projection) of a hexagon lie 
on two lines I, I' at the same z; the remaining sites are two 
pairs (lying on lines in the other diagonal direction). Hence, 
the loop product for that hexagon is r s = Yii£ S Vi — v l v l'- 
If we choose {v?} to alternate between adjacent lines in the 
same z-layer, we ensure that r r = — 1 for all hexagons as e.g. 
in Fig.^c). [Configuration (d) in Table H]is a layered state, 
with signs v% alternating in each layer of spins along [110] 
lines but all signs vg = +1 for the lines of spins in the [110] 
direction and hence not a harmonic ground state.] There is 
still an arbitrary sign choice in each z-layer, so this construc- 
tion gives 2 Z distinct ground states, where Z is the number of 
distinct layers (separated by Az = a/4 in terms of the cubic 
lattice constant a). This gives a lower bound of 0(L) for the 
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FIG. 1: [COLOR ONLINE]. Spin configurations on the pyrochlore 
lattice. Sites shown lie in a slice of thickness one unit cell, normal 
to the page, and all these configurations repeat in that direction. The 
hexagon and puckered octagon loops are superposed (dashed). The 
diagonal lines lie in the plane of the page, forming four families at 
different heights z (marked). Lines running in the x or y directions 
are tilted 45° from the vertical, and are shown broken where they 
connect to an adjacent slice. The lattice unit cell (outlined in (a)) 
contains 4 spins and has lattice constant a. Configurations (a) and 
(b) are "gauge" equivalent and have the highest possible harmonic 
zero-point energy T . Configuration (c), with its 32-spin magnetic 
cell outlined, is the simplest one of many that have the lowest jF h . 
[Energies for these configurations are in TableQ] 



harmonic-ground-state entropy (in a cube of volume L 3 ), so 
the entropy per spin is zero in the thermodynamic limit. 

We do not yet understand the full set of ground states; note 
that "gauge" equivalent states exist that are not given by the 
layer construction llq|. and the total entropy was proven to 
have an upper bound 1 16] of 0(Z In Z), only logarithmically 
larger than the lower bound. 

Discussion — In summary, I have computed the O(JS) 
spin energy term which (partly) breaks a classical degener- 
acy in the large-S" pyrochlore antiferromagnet with purely 
isotropic interactions. Various complex orderings are found 
in real pyrochlore antiferromagnets, most often explained by 
elastic distortions or dipolar interactions l2fill or - less inter- 
estingly - by nonnegligible second-neighbor exchange. It was 
nevertheless valuable to isolate the role of quantum fluctua- 
tions here, since materials may be found in which those per- 
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turbations are very small 12111 . or frustrated Heisenberg mod- 
els may be cleanly realized by cold gases in optical traps 1 22 ] . 

Two tricks - the Moessner-Chalker equation of motion, 
Eq. (|3 and writing the zero-point energy as the trace of 
a matrix, Eq. (0 - enabled an (uncontrolled) expansion of 
the O(JS) term of the energy to give the effective Hamilto- 
nian 7F h . The result l|6} is a sum of products of Ising spins 
around loops. The problem has an exact "gauge" degener- 
acy implying a ground state entropy, to harmonic order, of 
(at least) O(L). To fully resolve the degeneracy, as in the 
kagome case JTolfl5ll . anharmonic spin-wave calculations are 
in progress 1 23 ] . 

It was plausible that the quantum system shows greater or- 
der than the classical one; what is surprising is that the py- 
rochlore is less ordered than the kagome Heisenberg antifer- 
romagnet in the classical T — > case |8], but more ordered 
in the quantum large-S* case, since harmonic fluctuations re- 
move an extensive degeneracy in the pyrochlore but not the 
kagome case. 

The loop effective Hamiltonian has already been applied 
to other systems 12411 . A similar approach made it possible 



to determine the ground state within the large-TV, large-/? ap- 
proximation of the pyrochlore 1 25 ] . A more complete analytic 
derivation 1 16] of (^} matches numerical fits to ~ 1%; it also 
highlights the modes (fewer than in the kagome case) which 
have divergent fluctuations in the harmonic approximation. 

An effective Hamiltonian such as (|6} has value beyond the 
possibility (as here) that it leads us to an unexpected ground 
state. It can also be plugged in to define a Boltzmann distri- 
bution, such as exp(— /3jF h ), which at low but nonzero tem- 
peratures is more valid than the classical spin ensemble fl2ll . 
It also a gives basis on which to build more complete or more 
realistic models, by the addition of anisotropies, quantum tun- 
neling [26], or dilution I12u27ll . Notice how the collinear se- 
lection has provided a different route than "spin ice" to realize 
an effective Ising model in a pyrochlore system; furthermore, 
tunneling between collinear states [26] might realize ring ex- 
changes II28I1 (and the consequent U(l) spin liquid). 
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